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Abstract Q value, the ratio of charmless hadronic decays between 7’ and J/1, is evaluated to be 
(26.0 + 3.5)%. The dependence in evaluation is taken into account carefully, and several approaches are 
adopted to find out the correlation coefficient or covariance between different measurements. 
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1 Introduction 


From the perturbative QCD (pQCD), it is expected that both J/w and y' decaying into light hadrons 
are dominated by the annihilation of cc into three gluons, with widths proportional to the square of the 
wave function at the origin [1]. This yields the pQCD “12% rule”, that is 


Byrn _ Byrete- 


D. = = (12.3 + 0.7)% . (1) 


Breck. -DPBilsche 


The violation of the above rule was first observed in px and K** K^ + c.c. modes by Mark II [2], 
since then BES has measured many two-body decay modes of %’ which violate this rule [3]. There are 
some modes which are suppressed in 7)’ decays relative to 12% rule, like Vecter-Pseudoscalar (VP) and 
Vecter-Tensor (VT) decay modes; while there are others which are enhanced, like KLK? mode. There 
have been many theoretical efforts trying to explain the puzzle [4], however, none explains all the existing 
experimental data satisfactorily and naturally. 

Some phenomenological studies indicate [5, 6, 7] that the S- and D-wave charmonia mixing scenario 
provides us an unified model to explain both the suppressed and the enhanced decays of Y’. If the scenario 
of S- and D-wave mixing is correct, it should give correct relation between the branching fraction of J/~, 
yw’ and vj" decays to the same mode. Therefore by virtue of measurements of J/1» and 7)’, we can estimate 
the corresponding decay width of V)" to the same final states. 

A recent research [7] tries to utilize the available information from J/y and «)' decays to evaluate 
the corresponding charmless decay width of yw". For such study, it is crucial to know the proportion 
and extent of the suppressed and enhanced ' decays. A major impediment to doing estimation in a 
systematic manner is the dearth of 7’ branching fraction measurements. So an alternative approach, 
whose main idea is to obtain charmless decay proportion by subtracting the charmed decays, is employed 
to estimate Qa. In the light of the approach, a value Q, equivalent to Qa, which is defined as 


_ B(' > ggg + 799) (2) 
B(JI/b > ggg + 799) 


has been estimated based on experimental results to be Q} = (2347) % in Ref. [8] and Q7 = (24.0+5.6) 
% in Ref. [9]. 

It is clear that such Qg value is considerably enhanced with respect to Qp, let alone that of suppressed 
decay modes. Such estimation implicates that while some modes are suppressed in v/' decays, a substantial 
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fraction of 7’ decay modes is enhanced relative to J/wv, with the average enhancement factor of more 
than 2. 

It is also worthy of noticing that the difference between two previous Q,’s. Although Qi and Q? are 
estimated on the basis of the same values from PDG2002 [10], their central value differ by 5%, moreover, 
the relative error of these two estimations differ from each other by about 30%. Such large discrepancy 
should be understood. 

As Qg is important in understanding of charmonium decay, and the discrepancy between previous 
estimations is rather prominent, it is necessary to evaluate Q, again based on the new PDG [11]. 


2 Evaluation of Q value 


The estimation of Q, is based on an assumption that the decays of J/w and w” in the lowest order of 
QCD are classified into hadronic decays via three gluons (ggg), electromagnetic decays via virtue photon 
(y*), radiative decays into light hadrons via two gluons (ygg), and decays to lower mass charmonium 
states (ccX) [9, 12]. Thus, using the relation B(ggg) + B(ygg) + B(y*) + B(X) = 1, one can derive 
B(ggg) + B(ygg) by subtracting B(y*) and B(ccX) from unity. 

The calculated values of B(y*) and B(ccX), together with the values used to calculate them are 
summarized in Table 1. As regards to w', two faint branching fractions of final states y7(25) and 
he(1 P1) + X are neglected in our evaluation. Herein special attention should be paid to the dependence 
between various branching fractions. In fact, many approaches have been used to acquire the correlation 
coefficient or covariance between different measurements. 


Table 1: Experimental data on branching fractions for J/4» and w’ decays through virtual photon and 
to lower mass charmonium states used in our analysis. Most of data are taken from PDG [11], except 
for B(J/p, p! — 7* — hadrons), which calculated by product R - B(J/,' > wtp), where R = 
2.28+0.04 [13]. In estimating the errors of the sums, the correlations between the channels are considered 
meticulously. 


Channel B(Jfv) Bp’) ord. of v 
y* — hadrons — (13.44:0.33)96 (1.660.18)% Th) 
e*e- (5.93£0.10)%  (7.55+0.31) x 10-8 (D 
"" (5.88+0.10)% — (7.3 40.8) x 10? & 
TUTO (2.8 50.7) x 10? 9 
yo xX (25.2240.43)% GB02% 

Me (L3 £0.4)%  —(2.8:0.6)x10-7 O 
ata J[ (31.7£1.1)% (D 
1989 J [jb (18.8+1.2)% 9) 
nI (3.16+0.22)% @ 
n9 J fi (9.6::2.1) x 104 O 
YXco (8.640.7)% ® 
YXe1 (8.440.8)% © 
YXe2 (6.4+0.6)% © 
aX (13204) (77.4£2.5)% 


2.1 Estimation involving J/y decay 


In order to evaluate the charmless hadronic decay of J/1), we need the branching ratios Be, B,, B.» and 
B-n., and the correlation information involving these channels. Since the experiment for yne is unrelated 
with other experiments [11], so there is no correlation between Byn. and other branching ratios. As to 
other correlations, we will consider them one by one. 

The leptonic branching ratios of J/«» decay listed in Table 1 are actually the synthetic results from 
four experiment groups, whose measurements, according to PDG, are presented in Table 2. These results 


can classified into two categories: one from energy scan experiment, the other from the analysis of 
p > nT J/v, JI > €t. 


Table 2: Experimental data on branching fractions for J/w and v decays to lepton pair. Correlation 
coefficients between measurement of e*e^ and jut states are also listed. 


ord. Experiment Collab. B. (96) B, (96) T/T, Coeff. ( p) 
1 e'e scan BES [14] 6.09 + 0.33 6.08 + 0.33 1.00 + 0.07 0.167 
2 e'e scan MARKI [15] 6.9 + 0.9 6.9 + 0.9 1.00 + 0.05 0.927 
3 ain J| BES [16] 5.90 + 0.05 + 0.10 5.84 + 0.06 + 0.10 0.764 
4 ain J| MARKIII [17] 5.92 + 0.15 + 0.20 5.90 + 0.15 + 0.19 0.629 
5 ain J| CLEO [18] 5.945 + 0.067 + 0.042 5.960+0.065+0.050 0.997 0.014 0.534 


As to the scan experiment, let’s consider the relation 


Ra e. 3 
H Fa B. (3) 


By the virtue of error propagation formula, we have 


2. ,2 2 " 
Vg, Z Ve E Vy — 2peu Ve Vp , 


or 2 3 2 
V. tu, — VR 


; (4) 


Pép = QWeVy 

where vg,, , ve, vj, indicate the relative error for Rey, Be, By respectively. Since vp,,,, ve, vj, are readily 
obtained from values in Table 2, the correlation coefficients can be calculated by Eq. (4), and are listed 
in the last column in Table 2. 

As to the second kind of results, which were presented with statistic and systematic uncertainties, we 
notice that the systematic uncertainty is common for ete~ and utu” final states analysis. By virtue of 
the theory purposed in Ref. [19], if two measurements x; and zj, have the common normalized uncertainty 
of, and their measurements are reported as follows 


go = axle we. 
q^ = el eee) coy, 


where c; and gj are the uncommon uncertainties, then the correlation coefficient between x; and v; can 
be calculated by 
2 
ETTEY 
[e 
A umque (5) 
e; + (ziot)? Ajo; + (ajax)? 


In our case, c; and c; are statistic uncertainty and oy is 1.7% for BES [16] and 3.3% for MARKII [17]. 

In the following analysis, we directly use e? and p? (short for o7 and o7) to represent the covariance 
of branching ratio B, and B,. We know the reported results in PDG are the weighted average of several 
experiments, viz. 


Be = o aiBei , B, = 5 biB pi , (6) 
i=! i-1l 
where ; 
22 
€ _9 1 
ay = e2 , = 1/ 5 e , 
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with i(= 1,2,3,4,5) denotes the four groups of experiment measurements in Table 2. Then the transfor- 
mation relation between Be (B,) and Bei (Bpi) is expressed as 


Popes e 


BT-(Ba Boa Bes Ba Bes Bua Bu Bua Bus Das), 
here superscript 7' denotes the transposition of vector or matrix and 
sp a2 az as œ 0 0 0 0 .) 


with 


0 0 0 0 0 by bo b3 ba bs 


Using the relation in Eq. (8), the covariance matrix V for B, and B, can be formulated as [20] 


V-SVS!, (9) 

where V is the covariance matrix for Be; and Bpi, which has the following form 

e p1€1pA 
e p2€2 H2 
e3 paea H3 
2 
e e 
= 4 i Paap . (10) 
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Here naught elements have been suppressed and pj indicates four correlated coefficients listed in Table 2 
Notice the definition of a, b in Eq. (7), we obtain 
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i—1 
together with the definition of B, and B, in Eq. (6), we work out the synthetic results 


B. = (5.9392: 0.061)96 , 
B, = (5.9312: 0.064)96 , 


which are just the values given by PDG (refer to Table 1). At the same time, we also acquire the 
correlation between Be and B, as reflecting by V in Eq. (11). 

Next we consider the correlation between B.» and other branching ratios. As we mentioned in Table 1, 
we calculate the B,» by the product 

By» =R-By , (12) 
therefore B,» is correlated with 6, and with it alone. In order to figure out the correlation coefficient, 
we consider the relation 
By» 


B, ^C 
which is similar to that of Eq. (3), so the corresponding Eq. (5) can be utilized to work out the coefficient 
puy» between B, and B+, which is 0.524 for J/w decays and 0.987 for 7’ decays. 

Summarizing the fore analyses, we obtain the covariance for J/w estimation 


R= 


(13) 


B. B, By» Ban. 
Be c2 [m 0 0 
By Pepe yu T; Pyy* OpO-y* 0 
Viy = B. 0 Puy* Op Oy 05, 0 i (14) 
Bane R 0 0 Orne 


With above covariance matrix, we work out B(J/w — ggg) + B(J/ — ygg) = (73.48 + 0.59) %. 


2.2 Estimation involving 7’ decay 


Since 2002, the treatment of the branching ratios of the v and Xco,1,2 has undergone an important 
restructuring. 

When measuring a branching ratio experimentally, it is not always possible to normalize the number 
of events observed in the corresponding decay mode to the total number of particles produced. Therefore, 
the experimenters sometimes report the number of observed decays with respect to another particle in 
the relevant decay chain. This is actually equivalent to measuring combinations of braching frations of 
several decay modes. To extract the branching ratio of a given decay mode, the collaborations use some 
previously reported measurements of the required branching ratios. However, the values are frequently 
taken from the Review of Particle Physics (RPP), which in turn uses the branching ratio reported by the 
experiment in the following edition, giving rise either to correlations or plain vicious circles. The way to 
avoid various dependencies and correlations is to extract the branching ratios through a fit that uses the 
truly measured combinations of branching fractions and partial widths. This fit, in fact, should involve 
decays from the four concerned particles, Y’, yeo, Xe1, and Ye2, and occasionally some combinations of 
branching ratios of more than one of them. This is what is done since the 2002 edition [10]. According 
to such a global fit of the new PDG [11], we can acquire the correlation coefficients for most of channels 
(channels from @ to (9)) listed in Table 1, only the correlation coefficients for channels @, Q, and Q 
have to be considered separately. 

Similar to the analysis of previous section, B(@) (i.e. By») is only correlated with process B, and its 
correlated coefficient has been given in the previous section. As to B(Q) and B((3), their correlations 
with other quantities are neglected for their comparatively small absolute values. 

So based on our analysis and with the information from PDG, we obtain the correlation coefficients 
for y’ estimation 
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Notice the symmetry of covariance matrix, so the symmetric elements are omitted in above matrix. The 
total covariance can be calculated by 


© © 
oy = 5 5 PijTiOj s (16) 
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where the values of p;; are presented in Eq. (15). With Tin, we work out B(Y' > ggg) + B(y' > ygg) = 
(19.14 + 2.54)%. 


2.3 Q value 


Using the results of previous sections, the ratio of branching fractions of vj' to J/w decays into hadrons 

is given by 

_ B(' > ggg +799) 
BU/) > ggg + 199) 

The central value of above estimation is consistent with the forementioned estimations Q} [8] and Q7 [9], 

but with higher accuracy. 


Q; = (26.0 + 3.5)% . (17) 


The relation between the decay rates of ggg and gg is readily calculated in pQCD to the first order 
as [21] 
T(J/$ — ggg) 5 os(m.) 
Using o; (m.) = 0.28, one can estimate the ratio to be 0.062. A similar relation can be deduced for the w’ 
decays. If we adopt 6.2% as esimation for ratio B(ygg)/B(ggg), we obtain B(J/v — ggg) ~ (69.240.6)% 
and B(w' > ggg) c (18.0 + 2.4)%, while the “26.0% ratio" stands well for either ggg or ygg. Although 
Qg is considerably enhanced relative to Qa, it coincides with the ratio for KK? decay between y’ and 
J/b, which according to recent results from BES [22, 23] is 


D(J/v > ygg) 16 a (1 i 2920) (18) 


The relative higher Q; implicates that a substantial fraction of )j' decay modes is enhanced relative 
to J/1), or some modes which exist in y’ decays, are absent in J/w decays. 

In order to understand the discrepancy of relative error between Qi and Q5, we also preform the 
corresponding calculation without consideration of any possible correlation. First, we work out B(J/4» > 
999) + B(J/ > ygg) = (13.48: 0.54)96 and B(Y' > ggg) + B(j > ygg) = (19.14 x: 2.05)96, respectively, 
whose relative error is 1096 and 2096 smaller than those of dependently estimated results. Then using 
Eq. (17), we obtain Q, = (26.0 + 2.8)%, comparing with Q; = (26.0 + 3.5)96, the difference between two 
relative error is up to 2596, which is just at the same level of the discrepancy for previous estimations!. 
Because there is no detailed exposition about Q value estimation, we have no idea about what kinds of 
and to what extent the correlation has been considered in previous estimation, and we can merely point 
out that the different treatment of correlation may be the reason leading to the discrepancy of Q value 
estimation between Refs. [9] and [8]. 


3 Discussion 


It is worth to mention another approach for estimating Qna [9], which is to use the data on branching 
fractions for hadronic decays in final states containing pions, kaons, and protons that have already 
been measured for both J/q and v. They are rtn, K^ K-, pp, v^v 5, ppr?, 2(n^m-), 3(rtr), 
2(n*n-)n9, 3(1*n-)n9, ntn K* K-, 1*1- pp. Using the PDG data compiled in Table 3, we have 


11 
M Bv > fi) = (10.87 + 0.74) % 


i=1 


and 
11 


Y B 2 fi) = (1.01 + 0.19)% . 


i=1 
It follows that 
Qs 


EBW > f)/ X BUje > f) Bs 


i= 


1 
= (9.30 + 1.82)% . 


We notice ddthat the Q, not only differ from Qe greatly but is suppressed with respect to Qp as well. 
A remark is in order here. We know that most of multihadron final states in fact include sums of 
several two-body intermediate states. One thus observes a mixed effect which may deviate noticeably 
from the expected value of Q, even if a few of the two-body intermediates are severely suppressed. For 
example, the decay ij" — ntr K^ K- can proceed through intermediate state K *(892)9 K3 (1430)? -- c.c., 


1 As forementioned, the correlations of B(@) and B(@) with other quantities are neglected. Otherwise, if we assume 
that one of, or all of these two branching ratios are completely correlated with other quantities, we can obtain the following 


results: 
B(®) correlated: = (19.14 + 2.654)% = (26.0 + 3.62)96; 
B(@) correlated: = (19.14 + 2.579)% = (26.0 + 3.51)96; 
bothB(@) andB(@) correlated: = (19.14 + 2.693)% = (26.0 + 3.67)% . 


Table 3: Branching fractions for v and J/i decays [11], and Q, values are also calculated. 


Modes Qn (76) 
0707 E T 54+ 35 
0.237 + 0.031 1.0 + 0.7 42 + 30 
BB D 2.12 + 0.10 2.07 + 0.31 9.8 1.5 
3 07 E 21.2 X: 1.011 0.8 + 0.5 0.53 + 0.34 
0- BB D 1.09 0.09 1.4 + 0.5 12.8 + 4.7 
multibody decay 4.0 1.0 4.51.0 11.3 + 3.8 
4.0 + 2.0 1.5 + 1.0 3.7 3.1 
296 35+ 16 12.1+6.1 
33.7 + 2.6 3048 8.9 3: 2.5 
7.2 X 2.3 16+4 22.2 + 9.0 


8.0 + 2.0 13.3 + 3.5 


+: The weighted average of two new measurements, one from BES [24] (2.10 + 0.12)% and BABAR [25] 
(2.18 + 0.19)96. 


whose Q = 2.91.3, is greatly suppressed [26]. So the Qs is not exact ratio of y’ to J/w inclusive hadronic 
decay rates, but represents on average the ratio of the exclusive decay channels, as measured to date. 

In fact, the current results concerning the v^ decays is rather limited, even sum all charmless channels 
presented in PDG2004 [11], the total branching fraction is less 296. Such condition prevents us from 
accurately testing the Qa value. However, the estimation of Q value provides us some clues concerning 
the exploration of charmonium decay dynamics. Since many suppressed channels have been found, 
especially those such as pr which is greatly suppressed in w' decay, and if the Qp is really represent the 
averaged value of inclusive hadronic decay, the estimation of Qg indicates that either lots of enhanced 
decays are not discovered, or some particular decays only present at w’, or both cases exist. Therefore, 
systematic study of 7’ decays are anxiously awaited. 


4 Summary 


Estimation of Q value, the ratio of charmless hadronic decays between J/41» and v/', is very important to 
understand and explore the dynamics of charmonium decay. In our paper, the dependence in evaluation 
is taken into account carefully, and several approaches are adopted to find out the correlation coefficient 
or covariance between different measurements. The final evaluation of Q value is (26.0 + 3.5)%, whose 
accuracy is much better than prevoius estimations. 
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